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A Laboratory Plasma Experiment for Studying Magnetic 
Dynamics of Accretion Discs and Jets 
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ABSTRACT 

This work describes a laboratory plasma experiment and initial results which should 
give insight into the magnetic dynamics of accretion discs and jets. A high-speed 
multiple-frame CCD camera reveals images of the formation and helical instability of 
a collimated plasma, similar to MHD models of disc jets, and also plasma detachment 
associated with spheromak formation, which may have relevance to disc winds and 
flares. The plasmas are produced by a planar magnetized coaxial gun. The resulting 
magnetic topology is dependent on the details of magnetic helicity injection, namely 
the force-free state eigenvalue a gun imposed by the coaxial gun. 
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1 INTRODUCTION 



The accretion disc occupies a leading role in astrophysics, 
figuring prominently in young stellar objects (YSO), bi- 
nary star systems, and active galactic nuclei (AGN). An 
unsolved mystery is the origin of highly collimated bipo- 
lar jets and episodic flares associated with accretion discs. 
It was proposed some time ago that magnetic field dynam- 
ics can supply the necessary jet formation and collimation 
mechanisms (Blandford 1976; Lovelace 1976). Magnetohy- 
drodynamic (MHD) simulations have shown that jets are a 
natural consequence of a rotating disc in the presence of a 
magnetic field (e.g. Shibata & Uchida 1985). Jet structure, 
such as "knots," and also episodic behavior are observed in 
the simulations (Ouyed & Pudritz 1997; Goodson, Bohm 
& Winglee 1999; Nakamura, Uchida & Hirose 2001). The 
details of these models are unlikely to be tested by observa- 
tions anytime in the near future. Therefore, data from lab- 
oratory experiments could be very useful for this purpose. 
It should be noted that astrophysical jets have been com- 
pared theoretically with plasma guns (Contopoulos 1995), 
and that plasma experimentalists have interpreted coaxial 
gun plasma flows in the context of astrophysical jet mor- 
phology (Caress 1996). 

This work describes a new plasma gun based labora- 
tory experiment and initial results which should give insight 
into the magnetic dynamics of accretion discs and jets. This 
experiment is the first to use a plasma gun which explic- 
itly simulates the geometry and topology of a magnetically- 
linked star-disc system by using a co-planar disc-annulus 
electrode setup. The experiment reveals (1) the formation 
and helical structure of a magnetically-driven collimated 
plasma, similar to proposed models of astrophysical jets, and 
(2) plasma detachment which may be relevant for disc winds 



and flares and field-line opening in disc coronae. The result- 
ing magnetic topology depends on the details of magnetic 
helicity injection, namely the force-free state eigenvalue a 
imposed at the boundary. The plasmas in this work satisfy 
the MHD criteria globally (S » 1, pi « L, and Va « c, 
where S is the Lundquist number, pi the ion gyro-radius, L 
the plasma scale length, and Va the Alfven speed). It should 
be noted that S ~ 10 3 in the experiment, similar to that in 
MHD numerical simulations. 

This paper is organized as follows. Section II gives a 
qualitative description of how the magnetic dynamics of 
an accretion disc system are modeled in the laboratory. 
Section III describes the experimental setup. Section IV 
presents initial experimental results, and Sec. V provides 
a discussion. Finally, a summary is given in Sec. VI. 



2 MODELING DISC MAGNETIC DYNAMICS 
IN THE LABORATORY 

For non-relativistic accretion disc dynamics, such as in 
YSO's and some binary star systems, global magnetic field 
evolution can be described using MHD. The ideal MHD 
Ohm's law, 



E + U x B = 



(1) 



implies that an electric field E will be induced due to the 
sheared toroidal rotation of the disc in the presence of a 
magnetic field B. The _Br component creates an electric po- 
tential drop V between the central object and the accreting 
disc, and therefore magnetic helicity K will be injected into 
the disc magnetosphere at the rate dK/At — 2ipV (e.g. Bel- 
lan 2000), where ip is the poloidal magnetic flux emanating 
from the central object. 
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In the laboratory, a planar coaxial disc-annulus elec- 
trode setup simulates the topology and approximate geom- 
etry of the astrophysical star-disc system. However, rather 
than rotating the annulus to induce an _Er as in a real ac- 
cretion disc, an electric potential Vg Un is applied across the 
electrodes instead. The background poloidal B, which cor- 
responds to the field of the central object, is generated with 
an external coil. Again through Eq. (jjj), the background B 
and the applied voltage induce an E x B toroidal rotation 
in the plasma, achieving a magnetic helicity injection rate 
AK/At = 2f/jgu n Vgu n , where V'gun is the poloidal flux inter- 
cepting the inner disc electrode. This setup is commonly 
used to create spheromak plasmas (Rosenbluth & Bussac 
1979; Jarboe 1994; Bellan 2000). 

Therefore, the most basic magnetic dynamics near a 
star-disc system, namely the winding up of poloidal field 
lines generated by the central object by disc rotation, can 
be simulated in the laboratory by applying a voltage across 
coaxial electrodes in the presence of a background poloidal 
field. Of course there are limitations to the laboratory setup. 
For example, Keplerian sheared rotation and the existence 
and placement of a star-disc co-rotation radius cannot be 
exactly replicated. An experiment cannot hope to capture 
all the nuances and range of behaviors in all accretion disc 
systems. However, it can provide experimental data on the 
most fundamental plasma configurations which can result 
from a simple magnetically driven system. 



3 EXPERIMENTAL SETUP 

Figure |l| shows a side-view schematic of the experimental 
setup. A planar coaxial gun is installed on one end of a 
large vacuum chamber evacuated to 1.5 x 10~ 7 Torr. The 
chamber is large compared to the plasma so as not to affect 
the plasma evolution. The gun setup, enlarged in Fig. || 
includes (1) an inner electrode consisting of a 20.3 cm diam- 
eter copper disc (attached to end of blue re-entrant port), 
(2) an outer electrode consisting of a 50.8 cm outer diameter 
copper annulus (in green), (3) an external solenoid (in red) 
to produce a poloidal bias magnetic flux ifr linking the inner 
and outer electrodes, and (4) gas lines to deliver fast puffs of 
neutral hydrogen gas to the desired position of breakdown 
(adjacent to the electrodes). The gap between disc and an- 
nulus is 0.635 cm; plasma breakdown does not occur in the 
gap because the pressure-distance product there is too small 
to satisfy the Paschen breakdown criteria (Bellan 2000). A 
cylindrical coordinate system (R, 6, Z) is utilized, with the 
origin at the center of the inner electrode and +Z defined 
to be away from the electrode (toward the left in Fig. [2J). 
There are eight gas injection holes on each electrode, dis- 
tributed uniformly along the 9 direction of the inner and 
outer electrodes. The bias field is characterized by the pa- 
rameter V'gun, which is the total initial bias flux intercepted 
by the inner electrode; sample contours of ip are shown in 
Fig- § 

The plasma formation sequence is as follows. After the 
bias field and gas puffs are introduced, negative polarity 
high-voltage is applied to the inner electrode by discharg- 
ing a 120 /jF capacitor bank through an ignitron; the outer 
electrode is maintained at vacuum chamber ground. In the 
presence of the high-voltage, the neutral hydrogen breaks 




Figure 1. Side-view schematic of vacuum chamber and planar 
coaxial gun setup. 




Figure 2. Side-view schematic of planar coaxial gun, including 
gas feeds, external magnetic field coil, poloidal flux ip contours, 
and symmetry axis. 

down. The optimum path for breakdown is along vacuum 
field lines linking the inner and outer electrodes as shown in 
Fig. The gun voltage V^un drives a current 7 gun between 
the electrodes, twisting up the purely poloidal vacuum field 
and producing a Be- In an accretion disc, it is the disc ro- 
tation which achieves the same effect. The initally discrete 
flux tubes expand and eventually merge into an axisymmet- 
ric configuration. Throughout this process, magnetic helicity 
is being injected from the gun into the plasma. 

Experimental results to date have centered around 
time-resolved global imaging of the plasma evolution. Im- 
ages shown here were taken using a Cooke Corp. HSFC-PRO 
multiple-frame charge-coupled device (CCD) camera, which 
takes up to eight images per plasma discharge. The camera 
view is through the 20.3 cm window (labeled in Fig. Q) such 
that the gun electrodes appear on the right hand side of each 
frame. A false-color table is applied to the images for ease 
of viewing. Typically, the plasma is observed in unfiltered 
visible light. However, by using filters, it has been verified 
that most of the light emission is from neutral hydrogen line 
transitions. The exposure time of each frame is 20 ns and the 
interframe time is set typically to 1.5 /xs, which is on the or- 
der of an Alfven transit time. Additionally, a triple Langmuir 
probe is utilized to measure localized values of n e , T c , and 
floating potential. A Rogowski coil encircling the ceramic 
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break (see Fig. g) measures /gun, and a high- voltage probe 
measures V gun at the mouth of the re-entrant port (shown in 
blue in Fig. ^). Typical experimental parameters are as fol- 
lows: V eun ^ 4-6 kV, /gun » 70-130 kA, Vgun ~ 0.5-2 mWb, 
B ~ 0.2-1 kG, T c ~ Ti « 5-20 eV, n c ~ 10 14 cm" 3 , and 
global /? = 2fi nk(T c + Ti)/B 2 ~ 0.02-0.1. 



4 INITIAL EXPERIMENTAL RESULTS 

Because the plasma is relatively low-/3, it is reasonable as a 
first approximation to consider it as nearly force-free with 
plasma currents nearly parallel to the magnetic field, i.e. 



V x B ~ aB. 



(2) 



Equation (|2j) implies that a is constant along field lines but 
not necessarily across them. Integrating Eq. (H) over the gun 
surface, it can be shown that the force-free state eigenvalue 
imposed by the gun is cv g un = /Uo/gun/V'gun. Experimentally, 
both /gu n and V'gun can be adjusted to achieve a range of 
a gun values. Depending on the peak value of a gun , several 
distinct plasma configurations are identified. 
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Figure 5. A plot of » gun versus column length L for different 
plasmas, showing good agreement with the Kruskal-Shafranov 
condition for a current-driven kink instability. 



4.1 Collimated plasma: an analog for disc jets 

For low values of a gun , the formation of a collimated plasma 
is observed, as shown in Fig. In this plasma, c\ gU n peaks at 
approximately 66 m _1 . In Fig. [](a), gas breakdown has just 
occurred along eight discrete paths, each of which follows 
vacuum magnetic field lines and terminates at gas injection 
holes. In frame (b), the discrete arcs have expanded and be- 
gun coalescing. Magnetic reconnection is expected to occur 
along the Z-axis as the discrete flux tubes coalesce; the par- 
ticularly intense light emission there may be due to higher 
plasma density resulting from compression. In frames (c)- 
(d), a central column forms and expands in the ^-direction. 
In frames (e)-(g), the central column extends in length and 
persists for several Alfven transit times. In frame (h), the 
central column begins to break up. 

The bright structures in the images are expected to cor- 
relate with magnetic topology and current flow; this was 
verified previously in a similar experiment using direct mea- 
surements of magnetic field (Yee & Bellan 2000). This cor- 
relation is consistent with (1) the light emission being pre- 
dominantly from neutral hydrogen atoms excited by current- 
carrying electrons and (2) the fact that the plasma is low-/3 
and expected to be in a nearly force-free state in which J 
is nearly parallel to B. It should be noted also that H-H + 
charge-exchange time is estimated to be very fast (<C 1 /us) 
in this experiment, and therefore light emission may be rep- 
resentative of plasma ion dynamics. Filamentary structures 
can be seen inside the column, suggesting complex field 
topologies within the globally collimated structure. 

The cross section of the central column is clearly ob- 
servable in Fig. |3| If there are no collimating forces, clearly 
the length of the column would not be as elongated as ob- 
served, nor would the cross section be so uniform along the 
Z-direction. This result verifies that the necessary magnetic 
structure and collimation for a disc jet can arise from mag- 
netic forces associated with helicity injection. 



4.2 Helical instability: an analog for jet structure 

Higher a gun results in plasmas with unstable central 
columns, as shown in Fig. kl in which Q gun peaks at ap- 
proximately 71 m -1 . In Fig. We), a highly nonlinear helical 
perturbation appears in the collimated plasma. By alter- 
ing the CCD camera timing, it is possible to determine the 
characteristic growth time of the instability, which is about 
1.7 /its, similar to the characteristic Alfven time. The insta- 
bility appears to be an ideal current-driven kink mode. 

The Kruskal-Shafranov condition for ideal kink stability 
in a cylindrical column is (e.g. Bateman 1978) 

g cdgc = 2-rraBz/LBg > 1, (3) 

where a is the column radius and L the column length. As- 
suming that both the initial V'gun and the instantaneous / gun 
are fully contained inside the plasma column, and using the 
definition of a gU n and i/> gun , it is straightforward to show 
that the condition for stability is 

Qguni < 47T. (4) 

Figure |H| plots a gun versus L for a collection of central 
columns taken from both different discharges and differ- 
ent times from within the same discharge. Stable (tri- 
angles), marginally unstable (squares), and kinked (dia- 
monds) columns are plotted in a g u n -Z< space, along with the 
Kruskal-Shafranov stability threshold (dashed line), showing 
good agreement between experiment and theory. This result 
indicates that plasma instabilities, in this case an ideal kink, 
can give rise to macroscopic structure within jets. 

4.3 Detached plasma: an analog for disc flares 

Still higher values of a gun result in detached plasmas, as 
shown in Fig. [], in which a gU n peaks at approximately 
129 m" 1 . In this discharge, the plasma appears to detach 
from the electrodes in frame (d) and then propagates along 
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Figure 3. Images of plasma evolution (shot 1210; peak a gun 66 m 1 ) in which a plasma column forms and persists 
for many Alfven transit times, illustrating the magnetic topology required for an astrophysical jet. 
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Figure 4. Images of plasma evolution (shot 1233; peak a gv 
current-driven kink, develops on the ideal MHD timescale, illustrating one possible source of jet internal structure. 



71 m 1 ) in which a helical instability, likely a 



the Z direction at a speed of approximately 6 x 10 4 m/s, 
a fraction of the estimated Va- It is likely that a sphero- 
mak configuration is formed here; this was verified previ- 
ously with direct measurement of B on an experiment with 
a non-planar source but similar helicity injection (Yee & Bel- 
lan 2000). This result supports the idea of field line recon- 
nection above accretion discs, which can lead to disc winds 
and also episodic high energy flaring. 



5 DISCUSSION 

The results above show three distinct plasma configura- 
tions having accretion disc characteristics. All three con- 
figurations result from the same plasma formation process, 
the only difference being the peak value of the parameter 
Qgun = Holgun/ipgun- Figure |^ illustrates this dependence by 
placing different plasmas into Ig Un -ip SU n parameter space, 
with detachment at larger a gun , attached columns at lower 
a gun , and kinked columns near a cr it ~ 60-70 m _1 . 

This a gun dependence has many implications. Most im- 
portantly, it suggests that the plasma configurations asso- 



ciated with accretion discs and jets are related to Taylor 
relaxation theory (Taylor 1986), which is a description of 
how a plasma evolves as magnetic energy is minimized sub- 
ject to the constraint of constant magnetic helicity. This 
process can be cast as a variational problem, and it can 
be shown that the resulting magnetic field configuration 
satisfies Eq. (Q) with uniform a. In the case of a driven 
plasma like an astrophysical jet, Va will be non-zero, and 
magnetic helicity will flow from regions of high to low a, 
which tends to "relax" the plasma toward uniform a (Bel- 
lan 2000). The evolution of many laboratory plasmas, in- 
cluding spheromaks, can be understood in terms of Taylor 
relaxation. For example, it has been shown that a threshold 
Qcrit at the source must be exceeded in order to form the 
closed-field configuration of a spheromak (e.g. Yee & Bellan 
2000). This property is not surprising since analytic solu- 
tions of Eq. ^) show that B transitions from sinh-like to 
sine functions above a critical a (Bellan 2000). 

Well-known models of astrophysical jets have consid- 
ered separately the roles played by poloidal (Blandford & 
Payne 1982) and toroidal fields (Contopoulos 1995) in jet 
formation. However, the present work shows that jet struc- 
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Figure 6. Images of plasma evolution (shot 1181; peak ct gun ~ 129 m 1 ) in which the plasma detaches from the 
electrodes, illustrating the possibility of field-line opening in disc coronae. 
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Figure 7. A plot of peak / gun versus Vgun f° r a collection of 
plasma shots, showing the dependence of the plasma evolution 
on the peak value « gun = /^o-fgun/V'guni with a cr it ~ 60-70 m — 1 . 



of disc jets and flares, respectively. The onset of a helical 
instability in the plasma column is shown to be consistent 
with the Kruskal-Shafranov condition for an ideal current- 
driven kink. The magnetic topology depends on the force- 
free state eigenvalue a gun imposed at the plasma gun. It 
is argued that Taylor relaxation provides a useful descrip- 
tion of the magnetic structures of accretion discs and jets. 
These results demonstrate experimentally that the concept 
of magnetically-driven jets is a viable one. More quantita- 
tive characterization of the observed plasmas is planned. 
This will include direct measurements of the magnetic field 
via insertable probes, ion flow velocity along the collimated 
plasma via Doppler spectroscopy, and plasma pressure pro- 
files via Langmuir probes. 
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ture should result when the correct value of a occurs at 
the disc boundary, resulting in a non-arbitrary combination 
of poloidal and toroidal fields in the disc corona. Likewise, 
deformation of the jet structure will occur under different 
choice of a. In the experiment, the observed instability is 
consistent with an ideal current-driven kink. This instabil- 
ity has been observed in numerical simulations of AGN radio 
jets (Nakamura, Uchida & Hirose 2001) and proposed as the 
responsible mechanism for wiggled structure in the jets. 



6 SUMMARY 

A plasma gun based laboratory experiment approximates 
the boundary conditions and topology of a star-disc system, 
and it is shown experimentally that magnetic helicity in- 
jection with these boundary conditions leads naturally to 
both collimated plasmas and detached plasmas, suggestive 
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